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1. Introduction 


This report presents a study of the basic equations governing orbital 
transfers using long orbiting tethers. A very simple approximation to the 
general transfer equation is derived for the case of short tethers and low 
eccentricity orbits. Numerical examples are calculated for the case of in- 
jection into a circular orbit from a platform in eccentric orbit, and injection 
into eccentric orbit from a platform in circular orbit. For the case of long 
tethers, a method is derived for reducing tether mass and increasing payload 
mass by tapering the tether to maintain constant stress per unit area of tether 
cross section. Formulas are presented for calculating the equilibrium orbital 
parameters taking into account the mass of the platform, tether, and payload. 

2. Orbital Center of a Tethered System 

In order for a tethered system to be in equilibrium in a circular orbit, 
the sum of the centripital and gravitational forces on the system must be zero. 
If the system is gravity gradient stabilized in a vertical configuration so 
that all parts of the system move with the same angular velocity, then the 
condition of equilibrium is that 

l"2 V2 

f GMdm f A ? 

J — -J 0 ) 

'"1 '•i 

Equation (1) can be solved for the angular velocity which is 



dm r 


( 2 ) 


A particle at a position r defined by 


GM 2 

— = r.o)" 


experiences no net acceleration. Solving for r we have 

c 


r 


c 




( 3 ) 


This value of r^. will be considered to be the orbital center of the system. 

The motion of this part of the system is the same as that of a free particle at 
that radius. This definition gives a slightly different answer than computing 
the center of mass of the system. Equation (1) is valid in general. However 

it is more convenient to evaluate the contributions to the integrals of discrete 
masses using the expressions 



and 

In an eccentric orbit, the position of the orbital center will vary somewhat, 

but this effect can probably be neglected.* It may also be possible to neglect 

the mass of the tether or payload in certain siutations. The usual methods of 

calculating the center of mass may also be adequate for short tethers. 

The transfer equations given in the next section use the orbit of the 

center of the system as a reference in computing the post release orbit of the 

payload or base launching vehicle. The release of a payload can be considered 

as a double launch since both the payload and the launching vehicle go into 

different orbits after the release. If the tether has significant mass and 

remains attached to the launch vehicle, the transfer equations apply to the 

orbital center of the launch vehicle plus tether. 

n — 

This approximation is valid only for small eccentricity. 



3 . General Transfer Equations 



Figure la 


Figure 1 b 


Figure 1 . 


Launch of an end mass from a tether system, 
part (a) the launch is to a higher orbit. In 
(bj the launch is to a lower orbit 


*n Figure 1 , Ci and C2 are the extrema of the orbit of the orbital center 
of a tethered system. The derivations do not depend on whether Cj is greater 
than, less than, or equal to C2. Similarly Si and S2 are the extrema of the 
orbit of the satellite at the end of the tether. The separation of the 
satellite from the center of mass at release is which may be either positive 
or negative depending on whether the satellite being considered is at the upper 
or lower end of the tether. That is. Si may be greater than Ci as shown in 
Figure la, or Si may be less than Ci (Fig. lb). The separation of the satellite 
from the orbit of the center of mass is z^. on the other side of the orbit. The 
sign of £2 must bt the same as the sign of £1. We wish to compute S2 as a 
function of Ci, C2, and Si. 


The semi -major axis of the orbital center is 

= (Cj + C 2 ) / 2 

The semi-major axis of the orbit of the satellite after release is 

Sg = (Si + S 2 ) / 2 

The velocity of zhe center of mass at release is 
The velocity of the satellite at release is 



Since v^ can also be written as 



we can equate the right sides of equations (4) and (5) to obtain the equation 



Solving for a^ we have 

2 I - . i-i 

Si ■ a^ Ci^ Cl a^ 


( 1 . 2^ + -liLl 

"s " ^Si ■ ci3 acCi3 


Substituting a^ = (Si + S 2 ) / 2 and solving for S 2 . we have 



5 . 








$2 = -Si + 


2ac Sj 

2a J - S,^ 


2a ^ - Si^ 


( 10 ) 


The apogee $2 goes to infinity if 2a ^ = Si^ or Pz/^f. ~ ~ 1-2599. If a 


c * c 

,/a = 2^/'"^ = 1.2599. If l 
c ‘ “ c c 

is that of a 200 km orbit, ii must be 1710 km to achieve escape velocity. 

If the orbit of the payload is circular we can set Si = S 2 ^ S in 
equation (6) for S 2 to obtain 

2 


-S + 2 (- 


2Si + _Si_\-l 

Ci^ aCi^ 


J- - I (?. 

2S " 2 ^S 




2^ 



1 


a,C,- 



s' ■ a^Ci 


Substituting = (Ci + Cz)/2 and solving for C 2 we have 

Cl' C, - (Cl + C2) 2C2 

S ^ (Cl + C 2)/2 Cl + C 2 

Ci‘^ + C2Ci^ 2C2S^ 

Cl** = C2(2S3 - Ci3) 

C 2 = Ci 7'(2S^ - Cr^) (11) 

Using equation (11), we can calculate the required prigee height to achieve an 
orbit of radius S using a tether of length from the payload to the orbital 
C'^nter. The distance from the launching vehicle to the orbital center will 
depend on tne >'sldtive masses of the payload, launch vehicle and tether. 


4. Injection Into Circular Orbit 


Figure 2 shows a plot calculated using equation (11) for injecting a pay- 
load into a circular orbit from the Shuttle in an eccentric orbit. Figure 2 

plots against the perigee height of the orbital center of the system before 
release. Equation (9), which is &2 ' is reasonably v/ell satisfied in the 
plot. For the last point at = 44 km, the value of iz is 500 - 200 = 300 km. 
This is within about of the value given by equation (9), tiamely, iz ~ 

7^1 = 7 X 44 = 308 km. 

5. Injection Into Eccentric Orbit 

Figure 3 shows a plot calculated using equation (10) for injecting a pay- 
load into eccentric orbit from the Shuttle in a circular orbit. The distance 
«.i of the payload from the orbital center at 200 km is plotted against the 
apogee height of the payload after release. The approximate formula iz » 7^) 
is satisfied at the beginning of the plot. The curve becomes quite non-linear 
toward the end of the plot which is terminated at approximately synchronous 
altitude. The payload could be circularized at synchronous height using an 
apogee kick motor. The mass of the tether system necessary to reach synchronous 
altitude is substantial, and the payload small. Tapered tethers must be used 
for the reasons described in the next two sections. With a small payload the 
Shuttle orbit would not be greatly altered by the release. 

6 . Constant Diameter Tethers 

In previous sections it has been noted that tether lengths of 1250 and 
1700 km would be required to launch payloads to synchronous altitude or escape 
velocity from Shuttle altitude. In this section, the tether stress per unit 
area of cross section is computed as a function of tether length with no payload 
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Figure 2, Plot of perigee height of the ci;nter of mass in km vs. ti 
in km for inserting a payload into a circular orbit of 
500 km altitude. The first column lists 4i and the second 
column is the periger height of the center of mass of the 
system before launch. 
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Figure 3. Plot of the apogee height vs. i\ in km for launching a 

pa.vload into an eccentric orbit with the center of mass of 
the system at 200 km before launch. 



Figure 4. A tether of constant cross section attached to a 
launching platform. The orbital center of the 
system is at a distance ri from the center of the 
earth. 

In Figure 4, a tether of constant cross section and density p is connected 
to a platform at distance r^ from the center of the earth. The orbital center 
of the system is at rj. We wish to calculate the stress in the tether necessary 
to support the weight of the tether itself. The force on each element of the 
tether is the difference between the centripital force and the gravitation 
force. The total force F at the point is given by 

F =y* (dm r u)2 - GM dm/r^) (12; 

ri 


Substituting dm = p A dr where A is the cross sectional area, and setting F = 
C A where C is the stress per unit area in the tether we have 


or 


Z h - Q iK J (ro)2 - GM/r2)d» 


r\ 


r2 


Z = pj' (ruj2 . GM/r2)dr 


ri 


Performing the integration gives 

C = p 


i (r2‘ - r^^) + gM {- — 

t T2 ri _ 


(13) 
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Figure 5. Stress per unit area vs. (r? - rj) for a constant diameter tether. 

The length in km is listed in the first column, and the stress per 
unit area in dynes/cm^ at the point ri is listed in the second column. 


Figure 5 shows a plot of C In dynes/cm^ vs. (ra-ri) with ri = 200 km + 
where Rg is the radius of the earth. The angular velocity u is set equal to 

since ri is the orbital center of the system. To be in equilibrium the 
tether has to extend below r^ and be connected to a heavy mass to counterbalance 
the force on the part of the tether above ri. A value of 1.45 g/cc was assumed 
for p. The break strength of kevlar is about 2.7 x 10^° dynes/cm^. We see from 
the plot that tether lengths of 1250 and 1700 km are impossible with this material 
even without a payload. Depending on the safety factor used, tether length 
appear to be limited to several hundred kilometers. The next section discusses 
a way around this problem. 

7. Tapered Tethers 

In very long tethers considerable tether mass can be saved, and payload 
mass increased by tapering the tether so as to maintain constant stress per unit 
area along the tether. 



Figure 6. Tether cross section vs. position along the tether 
for maintaining constant stress per unit area. 




In Figure 6, two masses m^ and m 2 are connected by a tether of variable 
cross section A. We wish to derive an equation for the cross section A as a 
function of the distance r from the center of the earth such that the stress 
per unit area of cross section is constant. The tension is therefore 

T = c A 04) 

where c is the maximum safe stress for the tether material. The mass dm of a 
section of tether of length dr is 

dm = p A dr 05) 

where p is the density of the material. The gravitational force dF^ on the 


element is 


dF = -dm GM/r2 

g 


and the centrifugal force dF^ is 

dF = dm r 
c 


(16) 


(17) 


where is the orbital angular velocity. The net tension force dT on the wire 
element is 


dT = f dr 


(18) 


Differentiating equation (14) and substituting into (18) gives 


dT = c ^ dr 


(19) 


For the system to be in equilibrium, the sum of the gravitational, centrifugal 
and tension forces on each element of wire must be zero. Adding equations (16), 
(17), and (18) gives 


dA 

dF + dF + dT = -dmGM/r^ + dmrw^ + C dr = 0 
g c 


( 20 ) 


13 . 






Substituting equation (15) into (20), gives 

HA 

p A dr (-GM/r2 + roj^-) + c ^ dr = 0 , 

^ ^ (_GM/r2 + r o>2) dr (21) 

A C 


The differential equation (21) can be integrated to give A as a function of r. 
The result is 

jin A = - — (GM/r + i “^) + k (22) 

c ^ 

or 

A = t - ^ (GM/r + ^ r^ u^) + k (23) 

The integration constant k can be computed by specifying the cross section 
at some convenient point. For example if A = Aj at r = rj , then from equation 
( 22 ), 

k = iin Ai + -^ (GM/ri + ^ ri^ u^) (24) 

Specifying the cross section and therefore the tension [through equation (14)] 
at either end of the wire fixes the equilibrium angular velocity uo. For example 
the sum of the forces on mass 1 in equilibrium is 

-mi GM/ri^ + miriw^ + C Ai= 0 (25) 


Solving equation (25) for to gives 


4 


GM/rr 


c Ai 


(26) 


Note that c Ai must be less than the gravitational force for to to have a real 
solution. If the cross section ^2 is given at the upper end, the sum of forces 
on mass 2 in equilibrium is 


-m2 6M/r2^ + m2r2to^ - c A2 = 0 


(27) 


Solving for to gives 

0 ) = JgM/v2^ + 

1 ni2r2 

In this case, u 3 is always real, however it may be larger than the orbital 
angular velocity of an orbit above the earth's surface. 

The cross sectional area given by equation (23) has its maximum value at 
the "orbital center" of the system. That is. the value of r where the gravita- 
tional and centrifugal accelerations are equal. This can be shown analytically 
using equation (23). The maximum value of A occurs when the exponent in 
equation (23) has its maximum. The maximum is obtained by differentiating 
equation (22) with respect to r, which gives equation (21), and setting the 
result to zero. We have therefore 

-GM/r^2 -r r^ .2 = 0 

( 28 ) 

where r is the thickest part of the tether. The left side of equation (29) is 
m 

the gravitational acceleration, and the right is the centrifugal acceleration. 
Equation (29) is the equation for a free particle in a circular orbit of radius 

r . The net acceleration is away from the orbital center in both directions so 

m* 

that the system is supported from the orbital center, and the tether is thickest 
at this point. 

The total mass M of the tether obtained by integration of equation (2) 

is T2 

M =y* p A dr 

where A is given by equation (23). 

The net acceleration a 2 of the upper end of the tether is 


15 . 


&2 = - GM/ra^ + rz 


(31) 


The tension at the end is Ta = C Aa- This payload mass ma is therefore Ta/aa- 
Note that ra must be above the orbital center or aa and ma will be negative. 

A computer program has been written to evaluate the integral in equation 
(30) and plot the cross section as a function of r. Figure 7 shows a plot of 
tether cross section (cm^) vs. r (km) for the following parameters: ri = 6580 

km (202 km altitude), tether length = 1250 km, mi = 200 kg, p = 1.45 g/cc, 

C = 1.5 X 10^° dynes/cm^, and tether diameter = .06 cm at ri. The tension at 
the lower end is .424 x lO^ dynes. The orbital center is 600 km from the lower 
end. The tether mass is 776 kg. The tension at the upper end is .413 x lO^ 
dynes and ma = 214 kg. The orbital angular velocity is 1.0377204 x 10‘3 radians/ 
sec. The lower and upper end velocities are 6.8282 and 8.1253 km/sec. 

Figure 8 shows the tether cross section vs. r for a 1250 km tether attached 
to an 86 ton Shuttle at 202 km. The tether diameter at the Shuttle end is 3mm, and 
total tether mass is 6.9 tons. The payload mass is 160 kg. The tether is thickest 
at the orbital center which is 29.6 km from the Shuttle. The orbital angular 
velocity is 1.174905 x lO""* radians/sec. The lower and upper end velocities 

are 7.7309 and 9.1995 km/sec respectively. 

Figure 9 shows tether cross section for a 1700 km tether attached to the 
Shuttle at 202 km. The diameter at the Shuttle end is 4mm. The total tether 
mass is 13 tons and the payload is 32.6 kg. The orbital center is 53 km from 
the Shuttle. The orbital angular velocity is 1.168687 x 10"^ radians/sec. 

The lower and upper end velocities are 7.6900 and 9.6767 km/sec. 

Figure 10 has an 86 ton platform at the top of a 1700 km tether. The 
tether diameter is 4mm at the top. The total tether mass is 17 tons and the 
payload at the lower end is 187.5 kg. The ratio of payload to tether mass is 
about 4.4 times as great as in Figure 9. The orbital angular velocity is 
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Figure 7. Tether cross section vs. position along the tether for a 1250 km 
tether with 200 kg at the lower end and 214 kg at the upper end. 
The position in km starting at the lower end is listed in the 
first column, and the cross section in cm^ in the second column. 
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Figure 8. Tether cross section incm^ (second column) vs. position in km 
along the tether (first column) for a 1250 km tether attached 
to the Shuttle at the lower end. 
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Figure 9. Tether cross section in cm^ (second column) vs. position in km 
along the tether (first column) for a 1700 km tether attached 
to the Shuttle at the lower end. 
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Figure 10. Tether cross section in cm^ (second column) vs. position in 
km along the tether (first column) for a 1700 km tether 
attached to an 86 ton platform at the upper end. 


8.5360992 x 10’^ radians/sec. The lower and upper end velocities are 5.6618 
and 7.0679 km/sec. 

For any of the cases discussed above, heavier payloads may be launched 
by scaling all the masses and cross sections. As an alternative to increasing 
the Shuttle or platform mass, the Shuttle may be moved further from the orbital 
center so that less mass is required. 

G . Concl usions 

Equations have been derived for calculating the orbital parameters 
required for the Shuttle or Space Platform in order to launch payloads into a 
particular orbit using a tether of length n. The effect of the mass of the 
payload and tether have been considered and a method derived for tapering long 
tethers in order to reduce tether mass and increase payload mass. The general 
transfer equations developed make it possible to determine the post release 
orbit of both the upper and lower mass from the orbital parameter of the center 
of the system before release. This is important particularly in the case of 
launching heavy payloads from the Shuttle in low earth orbit in order to avoid 
re-entry, or assist in a planned re-entry. 


